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Abstract 

The QCD effects of twist-4 operators on the first moment of nucleon 
spin-dependent structure function gi(x,Q 2 ) are studied in the framework 
of operator product expansion and renormalization group method. We 
investigate the operator mixing through renormalization of the twist-4 op- 
erators including those proportional to the equation of motion by evalu- 
ating off-shell Green's functions in the usual covariant gauge as well as in 
the background gauge. Through this procedure we extract the one-loop 
anomalous dimension of the spin 1 and twist-4 operator which determines 
the logarithmic correction to the 1/Q 2 behavior of the contribution from 
the twist-4 operators to the first moment of g\(x, Q 2 ). 
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In the last several years there has been much interest in nucleon's spin structure 
functions gi(x, Q 2 ) and g2(x, Q 2 ), which can be measured by deep inelastic scattering of 
polarized leptons on polarized targets. Recent experiments on the nucleon spin struc- 
ture functions carried out at CERN |], || and SLAC [§, |J, have stimulated intensive 
theoretical studies on the nucleon spin structure functions [||. 

In the deep inelastic scattering, the perturbative QCD has been tested so far for 
the effects of the leading twist operators, namely twist-2 operators, for which the QCD 
parton picture holds. Now the spin structure functions would provide us with a good 
place to investigate higher-twist effects. Our purpose in this paper is to study the 
renormalization of higher-twist operators, especially the twist-4 operators, which are 
relevant for the first moment of gi(x,Q 2 ), that corresponds to the Bjorken and Ellis- 
Jaffe sum rules [^, |7]]. The anomalous dimension of the twist-4 operators determines 
the logarithmic correction to the 1/Q 2 behavior of the twist-4 operator's contribution 
to the first moment of g\. 

The first moment of the g±(x, Q 2 ) structure functions for proton and neutron turns 
out to be up to the power correction of order 1/Q 2 : 



where g\ (x,Q 2 ) is the spin structure function of the proton (neutron) and the plus 
(minus) sign is for proton (neutron). On the right-hand side, gA = Ga/Gv is the ratio 
of the axial-vector to vector coupling constants. Here we assume that the number of 
active flavors in the current Q 2 region is Nf = 3. Denoting {p, sl^j^iplp, s) = Aqs^, 
the flavor-5't/(3) octet and singlet part, a§ and ao = AS are given by 
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and A£ is related to the scale-dependent density AE(Q 2 ) which evolves as 

hence AS is the density at Q 2 = oo. Here we have suppressed the target mass effects, 
which can be taken into account by the Nachtmann moments f| . Note that taking the 
difference between and leads to the QCD Bjorken sum rule, the first order QCD 
correction of which was calculated in |9|, |K| [Tl| and the higher order corrections were 



given in 12, 13, 14, 15 



Now, the twist-4 operator gives rise to 0(1/Q 2 ) corrections |i~6 , |17| to the first 
moment of gi(x,Q 2 ). As can be seen from the dimensional counting, there is no 
contribution from the four-fermi type twist-4 operators to the first moment of gi(x, Q 2 ). 
The only relevant twist-4 operators are of the form bilinear in quark fields and linear 
in the gluon field strength. This is in contrast to the unpolarized case JT^|, where both 



types of twist-4 operators contribute. The common feature for the renormalization of 
higher-twist operator is that there appear a set of operators proportional to equations 



of motion, which we call EOM operators |L9, |20[ . And there exists the operator mixing 



among twist-4 operators which can be studied in the off-shell Green's functions where 
the EOM operators are inevitable. It should be emphasized that we have to keep the 
EOM operators to extract the physical observables like anomalous dimensions, which 
will be discussed later. 

The relevant operators in our case has the following properties; The dimension of 
the operators is 5 and the spin is 1. Its parity is odd and it has to satisfy the charge 
conjugation invariance. The flavor non-singlet operators are bilinear in fermion fields. 
Here we have to consider gauge variant EOM operators as well. 

Thus we have the following six operators which satisfy the above conditions: 



El = ^75 W - ^l 5 D a 0i/> - ^75 pD°i>, 

El = ^75 E-Y ^ + tf l5 E a + ^75 pD'i/j, (4) 
El = ^ 75 <9 CT fli) + ^75 #0^, El = ^757" # W + ^75 # iW, 

where = — igA a ^T a is the covariant derivative and G^ v = \e ^uapG 1 * 13 is the dual 
field strength. And we work with massless quarks for simplicity of the argument. 

Here one should note that not all of the above operators are independent, as in 
the case of twist-3 operators in g 2 (x,Q 2 ) plf , and they are subject to the following 
constraint: 

Rl = R° 2 + El (5) 

where we have used the identities, = ^{j^p} and [E^,E U ] = —igG^ u . Therefore 
any five operators out of (|j) are independent and they mix through renormalization. 

Here we take (R 2 , E±, E 2 , E 3 , E4) to be the basis of the independent operators. The 
only operator which actually contributes to the physical matrix element responsible 
for the Bjorken sum rule is R 2 . This twist-4 operator corresponds to the trace-part of 
twist-3 operator, (R r =3)a^ lf j, 2 — g^G cr { M1 7 Ai2 }?/> — traces, but there is no relation between 
the basis for the twist-4 and that for the twist-3 operators. 

We now study the renormalization of the operators. The composite operators, 0«, 
are renormalized by introducing the renormalization constants as 

(O i ) R = J2Z ll (0 J ) B , (6) 

j 

where the suffix R (B) denotes renormalized (bare) quantities. For the present basis 
we have the following renormalization mixing matrix: 
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The general features for the mixing matrix are the following [B^ [2(], P2| : (1) The 
counter terms for the EOM operators are supplied by the EOM operators themselves. 
This is because the on-shell matrix elements of the EOM operators ought to vanish. 
(2) A certain type of operators do not get renormalized. And if we take those operators 
as one of the independent base, the calculation becomes much simpler. (3) The gauge 
variant operators also contribute to the mixing. 

We compute Z^ by evaluating the off-shell Green's function of twist-4 composite 
operators keeping the EOM operators as independent operators. Thus we can avoid the 
subtle infrared divergence which may appear in the on-shell amplitude with massless 
particle in the external lines. Another advantage to study the off-shell Green's function 
is that we can keep the information on the operator mixing problem. And further, the 
calculation is much more straightforward than the one using the on-shell conditions. 

At the tree level, R2 operator contributes to the 3-point functions with quarks ip, 
ip and a gluon, in the external lines. So we consider the following one-particle 
irreducible (IPI) Green's function: 

Tt A = (0\T(OMp')A;m(p))\0) 1PI , (8) 

where the fields and the coupling constant involved represent the bare quantities. Here 
we employ the dimensional regularization (D = 4 — 2e) and take the minimal subtrac- 
tion scheme. The Green's functions are renormalized as follows: 



(r J fi = £z 2 ^3%(r >, (9) 

where Z2 and Z$ are wave function renormalization constants for quarks and gluon 
fields. We first present the evaluation in the usual covariant gauge. The one-loop 
radiative corrections arising from eight diagrams for R2 are represented as: 

(rf)l4oop = {l + ^[-^ 2 (^) + ^(G)]}(4f)tree 



1 g 2 1 



c 2 (G)](rlf) tree (io) 



where the quadratic Casimir operators are C*2(-R) = 4/3 and C^Cr) = 3 for QCD. In 
(|i~0|) , the tree-level Green's functions are given by 

(Itf )tree = iffWVn 

(rt A )tree = <?757.(p + p')^ a - ige^V^, 

(rg A )tree = -2<775<7<^+ /)^ a - 2<nb7,(p + p%T a 

+2757 CT (P + pV" - ige apa pl a ^T a 
(rlf)tree = -P757p(p + 

(Itf)tree = 91^+ ^ ~ <WY,(p + V%T a 

-gislaip + p') P T a + ige rjpa pl a 1 p T a (11) 

Here one can easily see that these five operators have their tree-level 3-point functions 
as linear combinations of four independent tensor structures. So in order to identify 
the counter terms properly as given in flT0|) we need to make use of the conditions for 
Zij extracted from the 2-point functions with x/j, ip in the external lines. 

Note that the tree- level tensor structure for R 2 , ige apa pl a ^T a ', appears also in those 
for Ei, E 2 and E±. Therefore, in order to extract the correct mixing-matrix element, 
it is crucial to keep the EOM operators. This feature is quite in contrast to the case 
of twist-2 operators, where we do not have to consider EOM operators at all. 

For the Green's functions of the EOM operators, we have additional Feynman 
diagrams due to the presence of the two-point vertices at the tree level. Further, 
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5 S 

the EOM operators like E% and E4 which are of the form E = ipB^, where B is 

dtp 

independent of fields, do not get renormalized: Z 44 = Z 55 = 1. 

To summarize we get the following result for the renormalization constants. (The 
detailed calculation will be discussed elsewhere PJ): 



Z ll — fC^-R); %12 — ^C 2 (R) — |C2(G), 

^13 — qG 2 (R) — g^M^), Z14 = ^C*2(G), 

z 15 — ~\d 2 (G), Z22 = \C2{R) + |C , 2(G), 

z 2 z = -\C 2 {R)-\C 2 (G), z 2i =\C 2 {G) (12) 

z 25 = \C 2 (G)i Z32 = — \C 2 (R) — ^C 2 (G), 

^33 = —2~C2{R) + gC2(G), Z34 = ~ ^^(G), 

^35 = — gC^i^?), ^44 = ^55 — 0. 

where we have written the renormalization constants as 



1 ,,- 

Zij = 5ij H — - Zij . (13) 
e I67H 



This result is in agreement with the general theorem on the renormalization mixing 
matrix discussed above. 

We now determine the anomalous dimension of R2 operator. In physical matrix 



elements, the EOM operators do not contribute [18] and we have 



(phys|(i^)fl|phys> = ^(physlORMphys) = (1 ~ z^^C 2 (R))(phys\(I%) R \phys). 

(14) 

Therefore the anomalous dimension j R2 turns out to be 

lR2 (g) = Z^-(Z^) = JL^ + O(g'), 7 ° 2 = 2z n = -fC 2 (R), (15) 
which coincides with the result obtained by Shuryak and Vainshtein |24| based on 



the background field method |25j in the coordinate space, where they discarded the 
contribution from the EOM operators by taking the on-shell quark external states using 
the equations of motion for massless quarks given by 

p^ = ^p = 0. 
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Here we also presents our result for the renomalization mixing of the twist-4 operators 
in the background field method We shall work with the momentum space. In this 
method we decompose the gauge field into classical background field and the quantum 
field as: 

Al = Af» + ^ 

and set up the Feynman rule, where we have an additional term in the three-gluon 
vertex contributing to this calculations. In the background field method, there 
appear only gauge invariant operators contributing the mixing through renormalization 



[ 27| . We take the independent operator basis to the three gauge invariant operators; 
R 2 , E\ and E 2 . Here we calculated the Green's function (§) with as the external 
gauge field. Taking into account the wave function renormalization constant of the 
background gauge field, we obtain the renormalization mixing matrix: 




1 + f C 2 (R)a/e \C 2 {R)a/e 



\C 2 {R)a/e 








1 + \C 2 {R)a/e -\C 2 {R)a/e 
-\C 2 {R)a/e 1 - \C 2 {R)a/e 




(16) 



where a = g 2 /167T 2 . This result leads to the same physically observable anomalous 
dimension of R 2 as given in (|15|) . 

Including the twist-4 effect the Bjorken sum rule becomes 




gl^Q^-g^Q 2 ) 
a s (Q 2 " 
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where ^3 is the reduced matrix element of R 2a , renormalized at Q^, which is defined 
for the general flavor indices, with t l being the flavor matrices, as 



R\ a = gi[>G av 'f1?il>, (p, s\R 2a \p, s) = f i8a 



0,..-,8). 



So far we have considered the flavor non-singlet part. Now we turn to the flavor 
singlet component. Here we note that there is only one non- vanishing independent 



gluon operator: G aa D^G^ a - This operator is equal to the flavor-singlet operator R° 2o . 
up to the gluon's equation of motion: 

G™D»G m = snh a G° a i/>. (19) 

So now we have only to take into account the mixing between = g^ r y a G' Ta ip and 

Eq = G au D»G m - gfty a & a i/,, (20) 

in addition to the previous results for the non-singlet part. The mixing between 
and Eg can be studied by computing the Green's function with two-gluon external 
lines, T^^ 4 , shown in Fig.l. Now we introduce the renormalization constant Z 16 as 

(R2) r - z n (R2) B + Z 12 (E 1 ) b + Z 13 (E 2 ) B + Z 1A (E 3 ) b + Z 15 (E 4 ) b + Z w (E g ) b . (21) 

From the diagrams of Fig.l, we get for the number of flavors Nf\ 

Zie = ~~ — ~ — o x -Nt, (22) 

hence we obtain the exponent for the singlet part 

7S 2N f 1 /32 2 A \ 

+ o N f ■ (23) 



2p 2p 3p A) V 9 3 
Including the twist-4 effects the first moment of cfc n (x, Q 2 ) becomes 



Jo 



(±L gA + l ag )(i _ ^1 + (a 2 )) + -AE(1 - 33 8Nf ^ + 0{a 2 s )) 







^[{±T^/3 + -/ 8 } — ^ +-/ -77^ J' ^ 



9Q 2 L L 12" 36 JSJ WW 9 JU V« S (Q 2 ), 
where /o, fz and /g are the twist-4 counter parts of do, 03 and as- /j's are scale 
dependent and here they are those at Q 2 ,. 

If we take into account the ghost terms in our QCD lagrangian, we get extra terms 
for the gluon EOM operator, which are expressed in terms of the ghost fields and satisfy 



the BRST invariance. In addition, there appears the so-called BRST exact operator 



T9|. EH, p2[ which participates in the operator mixing. However, it turns out that their 



contributions cancel with each other, and the final result does not change p3j . This can 



be more easily seen in the background gauge where we have only Eq for the additional 
independent operator and no ghost fields. 

Finally it should be noted that the matrix elements of the twist-4 operators f^s 
are considered to have ambiguities due to the renormalon singularity as discussed in 
the literatures [p8| . However, the exponents of logarithmic corrections to the 1/Q 2 
behavior, which we computed in the present paper, have definite values without any 
ambiguity. In case the Q 2 dependence of the moment fl24|) could be measured with 
enough accuracy in future experiments, we would be able to examine the presence of 
the twist-4 effects. 

We would like to thank K. Tanaka for valuable discussion. 
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Figure Caption 

Fig.l The Feynman diagrams for contributing to the mixing of R° 2 and Eq. 
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